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CARDINAL  INTERPOLATING  MULTIRESOLUTIONS  * 

ROBERT  MICHAEL  LEWIS  t 

Abstract.  The  basic  method  of  constructing  wavelets  is  by  means  of  a  multiresolution  approxi¬ 
mation  of  L^(R).  In  this  paper  we  present  a  class  of  multiresolution  approximations  for  which  the 
associated  scaling  function  has  a  simple  cardinal  interpolation  property.  We  present  the  construction 
of  such  multiresolutions  and  discuss  the  symmetry,  decay,  and  regularity  properties  of  the  associated 
scaling  functions  and  wavelets. 

Key  Words.  Wavelets,  interpolation. 

1.  Introduction.  Wavelets  have  recently  received  a  great  deal  of  attention  in  such 
areas  as  signal  processing,  image  processing,  and  the  numerical  computation  of  singular 
integral  operators.  In  its  simplest  form,  a  wavelet  is  a  function  ip  such  that  the  family 

ipjtk(x)  =  2^2ip(2^x  —  k)  j,k  €  Z 

forms  an  orthonormal  basis  for  L2(R).  Examples  of  such  ip  have  been  constructed  for 
which  the  ipj ik  also  form  bases  for  such  function  spaces  as  the  Sobolev  spaces  Hs( R). 
Moreover,  wavelets  naturally  exhibit  a  marked  degree  of  spatial  and  frequency  localiza¬ 
tion,  which  accounts  for  much  of  their  utility. 

A  simple  and  elegant  way  both  to  construct  wavelet  bases  and  to  compute  the 
wavelet  expansion  of  a  function  is  by  means  of  a  multiresolution  approximation  ([1],  [6], 
[8]).  A  multiresolution  approximation  M  of  L2(R)  is  a  sequence  (Vj),  j  6  Z  of  closed 
subspaces  of  L2(R)  such  that 

1.  Vj  c  Vj+1, 

2.  (J*L-oo  Vj  is  dense  in  L2(R)  and  Dj^-oo  Vj  =  0, 

3.  /(*)  6  Vj  <=►  f(2x)  €  Vj+1, 

4.  f(x)  €  Vj  4=^  f(x  —  2~^k)  £  Vj,  and 

5.  There  exists  a  function  g  £  V0  such  that  the  family  ( g(x  —  k))k£ z  forms  a  Riesz 
basis  for  V0. 

Of  paramount  importance  is  the  following  consequence  of  condition  (5): 

5'.  There  exists  a  function  (p  such  that  the  family  (<p(x  —  k))k^i  forms  an  orthonor¬ 
mal  basis  for  V0. 

The  function  (p  is  called  the  scaling  function  associated  with  the  multiresolution.  In 
this  notation  the  index  j  indicates  scale  while  the  index  k  indicates  spatial  location. 
The  scaling  and  translation  properties  (3)  and  (4)  imply  that  the  family 

<P,ik{x)  =  2j'2<P(2jx-k)  j,k£  Z 

forms  an  orthonormal  basis  for  Vj.  From  <p  we  construct  the  wavelet  ip. 

*  Research  sponsored  by  AFOSR-89-0363  and  the  State  of  Texas  under  contract  1059,  the  Texas 
Geophysical  Parallel  Computation  Project. 

*  Department  of  Mathematical  Sciences,  Rice  University,  Houston,  Texas,  77251-1892. 


1 


Interpolating  Multiresolutions 


2 


In  this  paper  we  will  study  multiresolution  approximations  of  L2(R)  for  which  the 
associated  scaling  functions  $  are  continuous  and  have  the  following  cardinal  interpo¬ 
lation  property:  For  k  6  Z, 


$(Jfc)  =  S(k)  = 


1,  if  k  =  0, 
0,  if  k  ^  0. 


We  will  call  a  multiresolution  for  which  the  scaling  function  $  is  continuous  and  satisfies 
this  interpolation  property  a  cardinal  interpolating  multiresolution.  Such  multiresolu¬ 
tions  are  also  discussed  in  [5]. 

The  main  result  of  this  paper  is  the  construction  of  examples  of  interpolating  mul¬ 
tiresolutions.  We  will  give  a  simple  condition  which  will  allow  us  to  construct  examples 
of  cardinal  interpolating  multiresolutions  and  we  will  study  some  of  their  properties. 
We  will  also  show  how  to  construct  cardinal  interpolating  multiresolutions  from  other 
multiresolutions.  We  will  construct  examples  of  cardinal  interpolating  multiresolutions 
from  the  spline  multiresolutions  of  Battle  and  Lemarie  [4]  and  from  the  compactly 
supported  multiresolutions  of  Daubechies  [2]. 

These  constructions  will  produce  families  of  cardinal  interpolating  multiresolutions 
for  which  the  scaling  functions  $  and  wavelets  'I'  are  real- valued,  exponentially  decaying 
functions  which  can  be  designed  to  have  an  arbitrarily  high  degree  of  smoothness  and 
an  arbitrarily  large  number  of  vanishing  moments.  As  we  will  discuss  in  Section  9, 
such  interpolating  multiresolutions  provide  an  attractive  and  natural  way  of  computing 
from  sampled  data  the  initial  projection  needed  for  Mallat’s  cascade  wavelet  transform 
algorithm  [7]. 

We  begin  in  Section  2  with  a  brief  review  of  some  properties  of  multiresolutions  that 
we  will  need.  In  Section  3  we  present  the  condition  which  allows  us  to  construct  cardinal 
interpolating  multiresolutions.  In  Sections  4 — 7  we  discuss  some  of  the  restrictions  that 
the  interpolation  property  places  on  the  multiresolution.  In  Section  8  we  show  how  the 
interpolation  condition  leads  to  some  interesting  relations  involving  the  values  of  the 
scaling  function  $  and  the  wavelet  tF  In  Section  9  we  discuss  the  vanishing  moments 
properties  of  interpolating  multiresolutions. 

The  main  results  of  this  paper  are  contained  in  Sections  11  -  13.  There  we  show  how 
an  arbitrary  multiresolution  can  be  used  to  construct  a  cardinal  interpolating  multires¬ 
olution  and  we  investigate  the  properties  of  the  cardinal  interpolating  multiresolutions 
so  constructed.  In  Sections  12  and  13  we  present  two  families  of  cardinal  interpolating 
multiresolutions. 

To  prevent  any  confusion  over  normalization,  the  definition  of  the  Fourier  transform 
we  use  in  this  paper  is 

A  roo 

fW  =  dxf(x)e—. 

J  —OO 


Its  inverse  is  then 
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For  a  generic  multiresolution  we  will  use  the  symbols  <f>  and  ip  for  the  associated  scaling 
function  and  wavelet.  We  will  reserve  $  and  $  for  the  scaling  function  and  wavelet 
associated  with  a  cardinal  interpolating  multiresolution. 

2.  Multiresolution  approximations.  In  this  section  we  will  review  those  fea¬ 
tures  of  multiresolution  approximations  (or  multiresolutions,  for  short)  which  we  will 
need.  For  a  fuller  discussion  of  this  subject,  we  refer  the  reader  to  [1],  [6],  and  [8]. 

Following  A.  Cohen  ([1]),  we  will  say  that  a  multiresolution  is  regular  if  the  scal¬ 
ing  function  (f>(x)  decays  so  quickly  that  (1  +  |a:|)n<?!>(:r)  G  L2(R)  for  all  n  >  0.  If  a 
multiresolution  is  regular,  then  its  properties  are  contained  in  what  we  will  call  the 
scale-transition  filter  H(u ).  This  is  a  function  which  describes  how  to  pass  between  two 
levels  Vj  and  VJ+1  of  the  multiresolution. 

Since  <f>( f )  G  V_!  C  V0,  we  have 

0^\o)  =  ^  hk<f>(x-k ), 

fc=— oo 

where 

(!)  hk  =  (|)  <f>(x  -  k). 

We  define  the  scale-transition  filter  H(io)  to  be 

OO 

■£  h 

k=  —  oo 

then 

(2) 

H( u)  is  C00  because  of  the  decay  rate  of  cf).  Since  —  k))^ z  forms  an  orthonormal 
family,  H(u)  satisfies 

(3)  H(  0)  =  1, 

(4)  l-ff(")|2  +  |R("  +  *)|2  =  1. 

As  a  consequence  of  (2)  we  have  the  fundamental  relation 

OO 

(5)  =  n  H{ 2-ku). 

k= 1 

From  (3)  and  (5)  we  see  that  <j>  has  unit  mean: 

(6)  i(0)  =  1. 

In  order  to  define  the  wavelet  ip,  we  first  define  the  filter  G(u>)  by 

(7)  G(  u)  =  e~iujH{u  +  7 r). 

The  wavelet  ip  is  then  defined  via 

w  ^“)=G(i)*(D- 
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3.  The  interpolation  condition.  We  define  an  interpolating  multiresolution  to 
be  a  regular  multiresolution  for  which  the  scaling  function  $  is  real- valued,  continuous, 
and  satisfies 

(9)  $(k)  =  6(k). 

We  should  remark  that  the  structure  of  multiresolutions  places  a  priori  restrictions  on 
the  kinds  of  interpolation  properties  we  might  consider.  One  property  of  multiresolu¬ 
tions  is  that  the  scaling  function  $  necessarily  satisfies 

OO 

y  —  k)  —  i. 

k=—oo 

(See  [8].)  In  particular,  we  must  have 

OO 

E  *(*)  =  !• 

k=— oo 

This  property  limits  the  interpolation  condition  we  might  hope  to  impose  on  the  scaling 
function.  In  view  of  this  restriction  the  condition  (9)  is  not  only  the  simplest  interpo¬ 
lation  property  but  the  most  natural  for  multiresolutions. 

We  will  derive  conditions  on  the  scale-transition  filter  H{uj)  under  which  the  associ¬ 
ated  scaling  function  4>  has  the  interpolation  property  (9).  We  begin  with  the  following 
well-known  formula  [8]:  Suppose  /,  /  €  L1,  and  let  F(u>)  —  ]C£L_00  f(u>  +  2kn).  Then 
F  €  L1[0, 27 r]  and 

OO  OO 

F(u>)  =  E  /(u.  +  2far)=  E  /(*)«“"• 

k=—oo  k=— oo 

In  the  case  of  an  interpolating  multiresolution  the  preceding  relation  reduces  to 

OO 

y  $(u>  +  2&7r)  =  i. 

k=  —  OO 

Now  suppose  that  H(u>)  is  the  scale-transition  filter  associated  with  the  multiresolution. 
Since  6(2a>)  = 

OO  oo 

E  *(2u,  +  2k7T)=  E  H(u  -f  kn)&(u:  +  kie)  =  1. 

k~ — oo  k=—o o 

Now  break  up  the  sum  over  even  and  odd  k.  Since  H  is  27r-periodic, 


H(u)  E  $(^>  +  2kft)  +  H(u>  +  7r)  E  +  tt  -f  2k-ir)  =  H(u)  +  H{u  +  tt). 


Thus  H  must  satisfy  the  relation 

(10)  H{u)  +  H(u>  +  7r)  =  1. 


Interpolating  Multiresolutions 


5 


The  condition  (10)  is  sufficient  to  generate  an  interpolating  multiresolution  under 
certain  other  assumptions  on  This  is  the  import  of  the  following  theorem,  which 

is  a  simple  elaboration  on  the  sufficient  conditions  concerning  the  construction  of  a 
multiresolution  from  a  scale-transition  filter  H(u)  in  [1]  and  [6]. 

THEOREM  3.1.  Suppose  H(u)  has  the  following  properties: 

(ii) 

H(uj)  ^  0  for  all  u>  €  [— tt/2,  tt/2], 

(12) 

H(d)  is  C°°  and  2n -periodic, 

(13) 

H(  0)  =  1, 

(14) 

\H(lo)\2  +  \H{lo  +  k)\2  =  1. 

Define 

OO 

$(w)  =  J]  H(2~ku). 

k=l 


A.  Then  $  defines  a  regular  multiresolution  of  L2(R).  The  function  $(a;)  is  the  Fourier 
transform  of  a  function  #(;r)  such  that  ($(x~k))k£ z  is  an  orthonormal  basis  for  a  closed 
subspace  V0  of  L2( R),  and  (1  +  |a:|)n<i>(;r)  €  L2(R)  for  all  n  >  0. 

B.  If,  in  addition,  $  is  integrable  and  H{u)  satisfies 

H{u>)  +  H(u  +  ir)  =  l, 


then  $  is  continuous  and  defines  a  cardinal  interpolating  multiresolution. 

We  will  need  the  following  lemma  for  the  proof  of  this  theorem.  This  lemma  is  a 
slightly  generalized  restatement  of  Lemma  2  in  [6]. 

LEMMA  3.2.  Suppose  m(cv)  is  a  continuous  2tt -periodic  function  satisfying 

m(ui)  -f  m{u  +  7r)  =  1. 


For  k  >  1  define  mk{u)  to  be 


mk(u>)  = 


m  (f)  ,  |  <  2fc7r, 

if  |o;|  >  2kir. 


Then  for  all  k  >  1  we  have 


/OO 

<komk(u)einw 

-OO 


27r,  if  n  =  0  , 
0,  ifn^O. 
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Part  A  of  Theorem  3.1  is  contained  in  [1]  and  [6j.  The  proof  of  part  B  uses  the 
techniques  from  the  proof  of  part  A. 

Proof.  We  must  show  that  <&(k)  =  6(k).  Since  <f>  G  L1,  $  is  continuous,  and  it 
makes  sense  to  speak  of  the  pointwise  values  §(k). 

Define 


3>fc(w)  = 


J 

r) H  (!) 

•••#(£),  ifM<2*7T 

1  0, 

if  |cu|  >  2fc7T 

Applying  Lemma  3.2  we  see  that 


/OO  A 

dco  ^(u^e* 

-OO 


27 v6(n). 


We  now  will  apply  the  dominated  convergence  theorem  to  conclude  the  same  about  $. 

The  conditions  (12),  (13),  and  (11)  imply  there  exists  a  lower  bound  c  >  0  such 
that  for  all  u  G  [— 7r,  tt], 

|$m|  > c- 

(See  [1]  or  [6]  for  a  proof.)  For  |o;|  <  2fc7 r  we  have 


=  **(«)*  (£) , 


so  for  |u>|  <  2fc7r, 
(15) 


$fc(w)  <  -  4>(u>) 
c 


Since  $k{ w)  =  0  for  |o?|  >  2fc7r,  (15)  holds  for  all  u>.  Since  we  assumed  that  $  is 
integrable,  we  may  apply  the  dominated  convergence  theorem  to  conclude  that 

$(n)  -  ~  r  dto  $(w)e,'nh'  =  6{n). 

Z7T  J- oo 


□ 


4.  Restrictions  on  the  real  and  imaginary  parts  of  H.  In  this  section  we 
will  summarize  the  conditions  we  need  to  impose  on  the  real  and  imaginary  parts  of 
the  filter  H(u )  in  order  to  produce  a  cardinal  interpolating  multiresolution.  These 
conditions  will  be  used  in  Section  11,  where  we  will  present  a  large  class  of  cardinal 
interpolating  multiresolutions. 

The  filter  H(ui)  must  satisfy  the  two  conditions 


(16) 


H(u)  T  H(lu  -f  7r)  =  1 


(17) 


|ff(u.)|2  +  |ff(k>-Mr)|2  =  l. 
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These  yield  \H(u)\2  —  R eH{uo)  =  0,  or 


(18) 


H{u) 


1 


1 

2 


Let  H(u>)  =  a{u>)  +  i  b(u>).  Then  (18)  leads  to  the  relation 

(19)  b2(to)  =  a(uj)  —  a2(u>). 

In  addition,  in  order  to  satisfy  (16),  a(o>)  and  b(to)  must  also  satisfy 


(20)  a(u>)  +  a(u)  -f  7r)  =  1 

(21)  b(u)  +  b(u  +  7t)  =  0. 

Finally,  if  $  is  real,  then  H(u>)  =  so  a  and  b  must  satisfy 


(22)  a(u>)  =  a(—  lo) 

(23)  b(co)  =  -b(-iv). 

Conversely,  it  is  easy  to  see  from  the  definition  of  4>  that  if  a  and  b  satisfy  the  two 
preceding  relations  then  $  will  be  real. 

5.  Missing  Fourier  coefficients.  The  scale-transition  filter  H{ui )  associated  with 
a  cardinal  interpolating  multiresolution  has  the  following  curious  property:  All  but  one 
of  the  “even”  Fourier  coefficients  of  H(u> )  vanish.  If  H(lo)  =  Y1T=- oo  bi-elkuj ,  then 


OO 

£  he'k“ 

k=—oo 


oo  oo 

+  £  hkeik^=  £  2 h,2ke'2k“ 


/c=— oo 


k=  —  oo 


=  1. 


Thus 


(24) 


b2k  = 


1/2,  if  A:  =  0, 
0,  if  k  ^  0. 


In  view  of  (1),  this  can  be  interpreted  in  terms  of  orthogonality  relations  between  cj)(x/2 ) 
and  the  family  ( cj)(x  —  k))kei. 


6.  Asymmetry  of  4>.  The  scaling  functions  associated  with  cardinal  interpolating 
multiresolutions  have  an  annoying  lack  of  symmetry.  In  fact,  the  interpolating  scaling 
function  $  cannot  be  symmetric  with  respect  to  the  origin.  Suppose  that  $(x)  is 
symmetric  with  respect  to  the  origin.  Then  the  Fourier  coefficients  hk  of  H(u>)  are 
symmetric  in  k:  hk  =  h_fc.  Accordingly,  H(u)  -  YlkL-oo  bk^ikuJ  is  purely  real,  so 


(lmH(u:))2  =  a(u>)  —  a2(u)  =  0. 


Since  a(ii>)  is  continuous  we  must  have  either  a  =  0  or  a  =  1.  Since  H( 0)  =  1,  it  follows 
that  a(u )  =  1,  which  contradicts  the  condition  a(u;)  +  a(o;  +  x)  =  1. 
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7.  Non-existence  of  compactly  supported  interpolating  multiresolutions. 


c 


=  !>• 


The  condition  (18)  restricts  the  values  of  H(u> )  to  lie  on  the  circle  {£  6  C  : 

As  a  consequence,  the  only  compactly  supported  multiresolution  with  the  interpolation 
property  $(k)  —  6(k)  corresponds  to  the  Haar  basis,  for  which  the  scaling  function  is 
not  continuous. 

To  prove  this,  we  first  note  the  following  fact:  If  p(£)  is  a  polynomial  and  |p(C)|  =  1 
for  all  |£|  =  1,  then  p(()  =  c(n  for  some  n  >  0  and  |c|  =  1  (see,  for  instance,  65, 
Part  VI,  of  [9]).  Now  suppose  that  H(u>)  has  a  finite  Fourier  expansion,  as  required 
by  Daubechies’  characterization  of  the  scale-transition  filters  associated  with  compactly 
supported  multiresolutions  [2].  Then  H(u>)  —  p(etLU)  for  some  polynomial  p ,  and  from 
(3),  (18),  the  aforementioned  fact  about  polynomials,  and  (24),  we  see  that 

-  \  = 


for  some  odd  n.  The  only  such  choice  of  n  that  leads  to  a  multiresolution  is 


H <">  =  ? + r‘"' 

which  generates  the  Haar  basis  (see  [1]). 

8.  Values  of  $  and  <5  at  dyadic  points.  For  cardinal  interpolating  multires¬ 
olutions  there  exist  simple  identities  for  the  values  of  $(A:2_J)  and  '&(k2~*)  for  any 
j,k€Z. 

Since  $(2w)  =  H( u>)6(u>),  for  j  >  1  we  have 


<l(2Jw)  =  H{  2j-1u;)$(2J'-1u;) 


Let  Qj(x)  =  <fr(2  *x).  For  j  >  1  we  have 


n  m  2*w) 


k~  o 


4>(u;). 


$j(w)  -  2>$(2iw)  =  2j 


n  m^) 


k= o 


Since  H  is  27r-periodic, 


OO 

y:  $j(o;  +  2kir)  =  2J 

k— — oo 


n  »(2^> 


k=z  0 


oo 

X  $(w  +  2kw). 

k=—oo 


Because  $  is  associated  with  a  cardinal  interpolating  multiresolution,  this  expres¬ 
sion  reduces  to  a  particularly  simple  form.  Since  $  satisfies  (9),  we  have 

oo  oo  j— 1 

X  <&{k2-j)e'k“  =  X  +  2^)  =  2J  n  H(2kL0). 

k=— oo  k=— oo  k— 0 
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From  this  relation  the  values  4>(fc2_J)  can  be  computed  in  terms  of  the  Fourier  coeffi¬ 
cients  of  a  finite  product.  This  identity  has  another  interpretation.  The  substitution 
u>  <—  k2~iio  yields 

oo  j— 1 

(25)  2~j  ^{k2~1)eik2~3^  =  J]  H(2~ku). 

A— — oo  A:— 1 


On  the  left-hand  side  we  have  a  rectangular  rule  approximation  to  <h(u;),  while  on  the 
other  side  we  have  a  truncated  product  approximation  to  $(u>). 

Similarly,  from  (8)  one  obtains  the  identity 

oo  j- 2 

J2  V{k2-j)eikw  =  VG(2i-lu)  H(2kco),  j>  2 

k=—oo  k= 0 


for  the  associated  wavelet  4/,  from  which  one  can  compute  the  values  4/ (A;  2  J)  and 
derive  a  formula  similar  to  (25). 

9.  Vanishing  moments.  In  Sections  12  and  13  we  will  show  how  to  construct 
cardinal  interpolating  multiresolutions  for  which  the  scaling  functions  $  can  have  ar¬ 
bitrarily  high  numbers  of  vanishing  moments.  In  this  section  we  relate  the  property  of 
vanishing  moments  to  the  scale-transition  filter  H(lo). 

For  an  arbitrary  multiresolution  the  number  of  vanishing  moments  of  the  scaling 
function  (j)  and  the  wavelet  ip  is  related  to  the  flatness  of  the  scale-transition  filter  H(u ) 
at  u  =  0.  First  note  that  if  H(lo)  <E  C°°,  then  <j>(u)  €  C°°  and  ip(uj)  e  C°°.  Moreover, 
from  (2)  and  (8)  we  see  that 

(26)  -^r-jr(O)  =  0  for  1  <  k  <  N  ^““^(0)  —  0  f°r  1  <  k  <  N,  and 

(27)  4r'(7r)  =  0  for  0  <  k  <  N  —  1  <=►  ^(0)  =  0  for  0  <  Jb  <  N  -  1. 

aw*  du>k 

Recall  also  that  H((2k  +  1  )7r )  =  0  for  all  k  6  Z  because  of  (3)  and  (4),  so  we  can  write 
H(u>)  in  the  form 


H(  u)  =  cosA  —m(u;), 


U) 


for  some  N  €  Z  and  some  function  m(u>). 

In  the  case  of  a  cardinal  interpolating  multiresolution,  H(u>)  also  satisfies  the  in¬ 
terpolation  condition  (16): 


H(u)  +  H(u>  -f  7r)  =  1. 
Consequently,  for  N  >  1  we  have 


drH 

duk 


(0)  =  0  for  1  <  k  <  N 


dkH 

dujk 


(7r)  =  0  for  1  <  k  <  N. 


(28) 
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This  leads  to  the  following  proposition. 

PROPOSITION  9.1.  Suppose  that  H(uj )  is  the  scale-transition  filter  for  a  regular 
cardinal  interpolating  multiresolution  I  and  suppose  that 


where  m(u>)  €  C°°. 
Then 


H(u>)  =  cosN  —  m(uj), 


u> 


J dx  xfc<f>(x)  =  0  for  1  <  k  <  N  —  1 
/ dx  x^^x)  =  0  for  0  <  k  <  N  —  1. 


Proof.  The  moment  integrals  make  sense  because  of  the  rapid  decay  of  $  and  'll. 
The  presence  of  the  factor  cosN(u/2)  means  that 


dkH 

du>k 


(7r)  —  0ior0<k<N  —  1. 


The  proposition  now  follows  from  (26),  (27)  and  (28).  □ 

One  useful  consequence  of  the  vanishing  moments  of  <f>  is  the  following  error  esti¬ 
mate.  Suppose  that  /  €  L2(R)  0  Cyv(R)  and  that  f^N\x)  is  bounded,  and  that  we 
have  a  cardinal  interpolating  multiresolution  for  which  the  scaling  function  $  is  expo¬ 
nentially  decaying  and  has  at  least  N  —  1  vanishing  moments  in  the  sense  of  Proposition 
9.1.  Given  Ax  >  0,  rescale  $  by  defining 


$A*(z)  =  (Ax)  5$  . 


The  family  —  kAx ),  k  €  Z,  is  an  orthonormal  basis  for  a  subspace  Vqx  of  L2(R); 

this  family  of  functions  can  be  used  to  define  a  multiresolution  of  L2  where  translation 
by  k  is  replaced  by  translation  by  kAx. 

Now  compute  the  coefficients  <4  of  the  projection  of  /  into  V£x,  the  subspace 
spanned  by  the  <&A;r(x  —  kAx),  k  €  Z.  Because  of  the  vanishing  moment  properties  of 
3>,  we  have 


<4  =  J  dx  /(x)$  Ax(z  —  kAx) 

r 

dx 


f(k Ax)  +  y;  fM{kAx)(x  -  kAx)"  + 


S‘ 

J  dx  [f(kAx)  +  i?(/V)(x)J  $Aa;(a:  —  kAx). 


$Ax(%  —  kAx) 


Recall  from  (6)  that  the  scaling  function  $  has  unit  mean.  Thus 

J  dx  f(k  Ax)$Ax(x  -  kAx)  =  (Ax)^  f(kAx). 
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Meanwhile,  for  the  remainder  term  RN  we  obtain 

|/ dx  R^N\x)$Ax(x  —  kAx)  <  K  J dx  |  —  kAx)\  |a;  —  . 

Since  $  is  assumed  to  be  exponentially  decaying,  say,  |$(a;)|  <  Ce~^x\  it  follows  that 

|/ dx  R^N\x)<f> Ax(x  —  kAx )  <  K'(Ax)N~*  . 

Comparing  the  interpolating  approximation  of  /  in  V^x, 

OO 

Y2  (Ax)5  f(kAx)$Ax(x  -  kAx), 

k=  —  OO 

and  the  exact  projection  of  /  into  the  subspace  V^, 

OO 

f(x)  die® AX{x  -  kAx), 

k— — oo 

we  see  that  the  first  approximation  interpolates  the  values  f(kAx),  while  the  error  in 
each  coefficient  is  bounded  by 

(Ax)*  f (kAx)  -  <4 1  =  O^Ax)^"*). 

Because  of  this  estimate,  interpolating  multiresolutions  provide  an  arguably  natural 
way  of  computing  from  sampled  data  the  initial  projection  needed  for  Mallat’s  cascade 
wavelet  transform  algorithm  [7].  We  can  interpolate  the  sampled  values  while  deviating 
from  the  exact  projection  with  a  controllable  error. 


10.  A  technical  lemma.  In  this  section  we  prove  a  lemma  which  we  shall  use  to 
prove  exponential  decay  for  the  multiresolutions  we  construct. 

Lemma  10.1.  Suppose  that  foruo  £  R;  9(lj)  is  2tt -periodic,  real- analytic,  |0(O)|  =  1, 
and  |0(u;)|  >  h  >  0. 

Then  0(C)  =  Id/tLi  $(2_fcC)  15  analytic  on  some  strip  V  =  {C  G  C  :  |  Im  (  |  <  r) 
containing  the  real  axis,  and  there  exists  a  constant  A  such  that  for  all  io  +  irj  €  T, 
|0(u;  +  iij)\  <  A  |0(w)|. 

Proof.  Since  9(u>)  is  real-analytic  and  27r-periodic,  9(f)  is  analytic  on  some  strip 
n  =  {(  €  C  :  |  Im  C  I  <  R}.  Because  |$(0)|  =  1?  if  K  CC  is  compact,  then  we  can 
find  Ck  >  0  such  that  for  all  (  €  K, 

|0(2-fcC)|  <1  +  CW2“*. 

This  insures  that  the  infinite  product  defining  0(C)  converges  uniformly  on  compact 
subsets  of  fi,  so  0  is  analytic  on  0. 

To  establish  the  bound  on  0(C),  we  first  define  for  u>  €  R 


A(w)  = 

A(w)  = 


9(u) 

OO 

E2-fcA(2-M- 

k=  1 
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Since  0(l>)  is  27r-periodic,  9' (to)  is  bounded  and  the  sum  defining  A(u>)  is  uniformly 
convergent. 

Since  A(tu)  is  27r-periodic,  we  may  write  A(cj)  =  YT=-oo  ^k^kw  ■  Because  of  the 
analyticity  of  9,  A(£)  is  analytic  on  fi,  so  for  some  a  €  (0, 1)  and  (3  >  1,  independent  of 
k ,  we  have  the  bound 

|Afc|  <  (3e~^. 

With  these  definitions, 

©v) = (£  n  = A(«)e<«). 

By  Leibniz’  rule,  for  n  >  1  we  have 

0(n)M  =  E  (  ”  ~k  1  )  A^)(u;)0^_1_^ (u). 

We  will  now  bound  0^nl(u;).  We  claim  that  for  some  constant  7  >  0,  independent  of  to, 

(29)  |e<”)H|<„!-2l|eM|. 

First  note  that 

CO  OO 

£  (ik)n \keik“  <  £  P\k\ne~aW. 

k=—oo  k=—c o 

From  the  identity 

f°°  xne~ax  —  n‘ 

Jo  ~  an+1 

and  Stirling’s  formula  we  can  derive  the  estimate 
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Then 


e(n)H  <  i;(n,1)|A<»||0(“-,-‘)M 


"  1  (n  -  1)!  7fc!  7"  1  k(n  -  1  -  k) I 
-  ^  kl(n  -  1  -  k)\  «*+*  an~k  1 

n— 1 


l — n  U. 


7 


k= 0 


which  establishes  (29). 
Then  for  any  u;  G  R, 


10(01 


A  0(n)M 


n=0 


n! 


/  oc 


<(E=rlC—  liieMI. 


\  z“'  an+1 

\n=0  u 


If  |C  —  a; |  <  <2/7,  then  the  latter  infinite  sum  is  convergent  and  bounded  independently 
of  to.  Thus  there  exists  a  strip  f  =  {(  G  C  :  |  Im  ( |  <  r}  containing  the  real  axis  and 
A  such  that 


|0(u;  +  it]) |  <  A  |0(u>) | 


for  all  o>  -f  ii)  G  R  □ 

11.  A  method  to  construct  interpolating  multiresolutions.  In  this  section 
we  give  a  method  of  constructing  an  interpolating  multiresolution  I  given  another  mul¬ 
tiresolution  M.  Recall  from  Section  4  that  for  an  interpolating  multiresolution  I  the 
filter  =  a(u>)  +  ib(u>)  is  determined  once  we  define  its  real  part  a(u>).  The  real 

part  a(o>)  must  satisfy 


a(uj)  +  a(u>  +  7r)  =  1. 

Happily,  examples  of  such  functions  a(u)  are  associated  with  any  multiresolution  M. 
Suppose  that  is  the  associated  scale-transition  filter;  then 

I#mMI2  +  |#m(^  +  tt)!2  =  1. 

Thus  we  might  take  as  a  candidate  for  a(cu)  the  function  |//m(‘*;)|2- 

Such  a  choice  of  a(u>)  has  the  following  attractive  feature.  Suppose  that  (f)  is  the 
scaling  function  associated  with  the  multiresolution  M.  Then  from  (18)  we  see  that 

\Hi(u)\  =  a^(u>)  =  |#mM|  • 


OO  OO 

®m|  =  n  |i/i(2-*^)i = n  |aM(2-‘w) 

k=1  k=  1 


(30) 

Consequently, 
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The  scaling  function  of  the  interpolating  multiresolution  I  is  as  smooth  as  the  scaling 
function  of  the  multiresolution  M  used  to  construct  it,  at  least  as  measured  by  the 
decay  of  the  Fourier  transform. 

In  the  same  sense,  the  wavelet  of  the  interpolating  multiresolution  I  is  as  smooth 
as  the  wavelet  of  the  multiresolution  M  used  to  construct  it.  The  wavelets  are  defined 
by  (8)  in  terms  of  the  filters 

GmM  =  e-'“HM(u  +  7r ) 

Gi(u)  =  e-luJH1(u>  +  tt). 

From  (30)  we  see  that  |£?i(u;)|  =  |£?m(^)|5  so 

As  we  shall  see,  interpolating  multiresolutions  so  constructed  inherit  most  of  the 
properties  of  the  multiresolution  used  to  construct  them.  This  gives  us  tremendous 
latitude  in  the  design  of  interpolating  multiresolutions.  We  will  present  two  classes  of 
interpolating  multiresolutions,  one  based  on  the  spline  multiresolutions  of  Lemarie  and 
Battle  and  one  based  on  the  compactly  supported  multiresolutions  of  Daubechies.  In 
both  cases  the  resulting  scaling  function  $  and  wavelet  $  are  exponentially  decaying 
and  have  the  same  smoothness  as  the  multiresolutions  used  in  their  construction. 

There  is  one  restriction  that  arises  in  the  construction  of  an  interpolating  mul¬ 
tiresolution  from  an  arbitrary  multiresolution.  This  restriction  arises  in  order  to  assure 
sufficient  smoothness  of  We  want  H\(u)  to  be  C°°  in  order  that  the  associ¬ 

ated  scaling  function  and  wavelet  are  rapidly  decaying  as  in  the  definition  of  a  regular 
multiresolution  in  Section  2.  In  the  examples  we  will  present,  Hy^{io)  has  the  form 

|#m(^)|2  =  cos27V  J  |m(u;)|2 

where  |m(u;)|2  is  an  even,  27r-periodic  function  and  |m(w)|2  >  0  for  all  u>  E  R.  Now 
construct  a(u>)  as  described  above: 

a(u)  =  \Hm(u)\2  =  cos2N  ^  |m(u>)|2  . 

From  (19)  and  (21)  we  have 

b2(uj)  =  a(u)  (1  -  a(w))  =  a(to)a(v  -f  t)  =  cos2/v  ^  sin2N  ~  |m(u>)|2  |m(u;  +  7r)|2  . 

Lj  Ld 

We  want  to  define  b(uj)  so  that  it  is  continuous  and  satisfies  b(u>)  +  b(u>  -f  7r)  =  0. 
Accordingly,  we  must  define  b(u>)  to  be 

sinA  (w/2)  cosA  (u;/2)  |m(cu)|  |?7r(tu  +  7r)| ,  N  odd 

a(u>)  sinw(u;/2)  cosn(lc/2)  jm(u;)j  j??i(  u  +  7r)  | ,  N  even, 


b(co)  = 


Interpolating  Multiresolutions 


15 


where  <r(u>)  is  the  27r-periodic  function 


a(u) 


1,  0  <  Ul  <  7T, 

—  1,  7T  <  U)  <  2r. 


If  N  is  even,  then  H(u)  cannot  be  C°°.  This  limits  the  decay  of  the  associated  scaling 
function  $  and  the  wavelet  we  cannot  expect  $  and  $  to  have  more  than  polynomial 
decay.  For  this  reason  we  will  limit  our  attention  in  the  remainder  of  this  paper  to  the 
case  where  N  is  odd.  The  constructions  we  present  can  be  carried  out  in  the  case  where 
N  is  even,  however,  yielding  multiresolutions  that  are  not  regular  in  the  sense  of  Section 
2.  Numerical  calculations  indicate  that  the  difference  in  the  asymptotic  rate  of  decay 
between  the  cases  N  even  and  N  odd  is  not  really  noticeable. 

THEOREM  11.1.  Suppose  that  f{x)  is  an  real-valued,  exponentially  decaying  scaling 
function  such  that  (1  +  |cj|)1+ct^  €  L°°(R)  for  some  a  >  0,  and  the  associated  scale- 
transition  filter  has  the  form  =  cosn(cj/2)  m(ui),  where  N  is  odd,  m(co)  is 

2ir -periodic,  |m(u)|  0  and  |m(a>)|  is  even. 

Then 

(31)  H\(uj)  =  (cosN  ~  |m(u>)|  +  i  smN  ^  |m(u>  +  tt)|^  Fm(w)  = 

defines  a  regular  multiresolution  with  scaling  function  $  and  wavelet  t[/  with  the  following 
properties: 

1.  Reality:  $  and  <P  are  real-valued. 

2.  Smoothness:  $  and  $  are  as  smooth  as  <f>  and  in  the  sense  that  for  u>  G  R. 

|$(u>)  = 

'F(uj)  =  . 

If  m  is  an  integer  for  which  0  <  m  <  a,  then  $  and  *P  are  Cm  functions. 

3.  Interpolation:  $  is  continuous,  and  <&(k)  =  6(k). 

4-  Exponential  decay:  $  and  *P  are  exponentially  decaying.  If  m  is  an  integer  for 
which  0  <  m.  <  a,  then  for  some  C  >  0  and  e  >  0, 

|4>(z)| ,  |4>'(x)| . . . . .  |*<”>(x)| ,  |4(x)| ,  !*'(*)[ , . . . ,  |*<”IP)|  <  Ce- 1*1. 

5.  Vanishing  moments:  $  and  'P  have  vanishing  moments: 

J dx  xk<f>(x)  =  0  for  1  <  k  <  N  —  1 

J dxxk'S>(x)  —  0  for  0  <  k  <  N  —  1. 


Proof.  First  we  must  check  that  defines  a  multiresolution. 

i.  Hfiuj)  is  2-7r-periodic  and  C°°  because  of  the  hypotheses  on  m(u). 

ii.  Since  //m(0)  =  1,  m(0)  =  1,  so  i/j(0)  =  1. 


Interpolating  Multiresolutions 


16 


iii.  Note  that 

|#(a;)|2  =  cos2/v  ^  |m(a;)|2  +  sin27V  —  |m(u>  +  7r)|2 

Li  Li 

-  +  \Hm{u  +  7r)|2  =  1. 

Thus 

|//lH|2  +  +  7T)j2  =  |i/M(^)|2  +  \Hm(u  +  ?r)|2  =  1- 

iv.  Since  \Hi(u>)\  =  cosN  j  |m(cd)|,  we  see  that  \Hi(u>)\  /  0  for  u  G  [— 7t/2,  7t2]. 
Theorem  3.1  then  tells  us  that  H\{u>)  defines  a  multiresolution. 

Now  we  proceed  to  the  proof  of  points  1 — 4. 

1.  To  prove  that  $  and  vf  are  real,  it  suffices  to  show  that  H\{io)  =  Hi(-ui).  Since 
(j>  is  real-valued,  it  follows  that  Hm(u)  =  #m(— whence  m(^)  =  ra(  —  a;).  Then 

Hi(-u)  =  .  .  . . (cos^  ^  |m(-u;)|  -  i  smN  ^  |  m(-u>  +  tv)\) 

m(—u>)  \  2  2  J 

By  hypothesis,  |m(u;)|  is  even  and  27r-periodic;  hence  |m(  —  u>  -f  7r)|  =  |m(u;  +  7r)|.  Thus 

-  Hi(uj). 

2.  We  have  shown  that  |LTi(u)|  =  so  part  2  immediately  follows. 

3.  Since  ^  €  L1,  part  2  tells  us  that  $  €  L1,  so  $  is  continuous.  To  establish  the 
interpolation  identity  $(k)  =  6(k)  we  need  to  check  that  H\(u>)  -f  H\(u  -f  n)  =  1.  We 
have 

-f  H\{u  -f  7r)  =  cos27V  ^  |m(u;)|2  +  isinN  ^cosN  —  |m(u;)|  |m(u>  +  7r)| 

Li  Li  Lu 

,  2/v  ^  i  /  \i2  .  ,  US  j\t  LO  .  /  \ii  /  \i 

+  sin  —  |m(u>  +  7r )  j  —  isin  —cos  —  \m{u>)\ \m[u>  +  7r)| 

Li  Li  Li 

=  \HM(u)\2  +  \HM(u;  +  ir)\2  =  l. 


4.  We  will  use  the  following  facts  to  prove  the  exponential  decay  (see  [3]). 

a.  If  |u(x)|  <  Ce~^x\  then  there  exists  a  strip  in  C  containing  the  real  axis  on 
which  u(Q  is  analytic. 

b.  If  u(C)  is  analytic  on  some  strip  T  in  C  containing  the  real  axis  and  if  u(()  = 
0(lC|-^1+/3*)i  P  >  0  as  £  — >  oo,  (  £  T,  then  for  some  C  and  e  >  0,  |u(x)|  < 
Ce~s  >*l. 

We  have 


$(w)  =  J]  6(2~ku)  cos2 


m( 2  kco)  = 


jj  9(2-ku)m(2~kLo) )  . 


Because  (f>  is  assumed  to  be  exponentially  decaying,  9(u>)  and  m(u )  are  real-analytic. 
Lemma  10.1  then  tells  us  that  0(C)  =  nr=i  ^(2  fcC)m(2  fcC)  is  analytic  in  some  strip 
T  =  {C  €  C  :  |  Im  C  |  <  r},  and  there  exists  A  such  that 


l«(u  +  i,)|  <  A  II  «(2-‘C)m(2-tC) 

k~l 


A  II  m(2-*C) 


k=  1 
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Meanwhile,  there  exists  B  such  that  for  all  (  =  uj  +  it]  E  T,  |cu|  >  1, 

-TV 


<  B 

\2/  2 

sin 


2TV 


UJ 


Consequently,  for  all  (  —  uj  +  ir/  €  T,  |u;|  >  1, 


$(C)|<  AB 


UJ 


-N 


uj 


-I  sm2N  — 

2)  2 


JJ  m( 2  ku ) 


k- rl 


AB 


4>(u) 


Thus  we  conclude  that  4>(()  is  analytic  on  the  strip  T  containing  the  real  axis  and 
<1(0  =  0(|C|-(1+a))  as  £  — >  oo,  so  $  is  an  exponentially  decaying  function. 

The  decay  of  consists  of  noting  that  for  0  <  k  <  m,  £fc$(£)  is  analytic 

on  T  and  is  0(|(|~*1+/3'),  0  >  0,  as  £  — >  oo  for  £  £  T.  The  proof  of  the  exponential 
decay  of  <ldm)  follows  the  same  lines. 

5.  This  is  simply  proposition  9.1.  □ 

From  formula  (31)  we  can  discern  something  of  the  relation  between  the  original 
scaling  function  <j>  and  the  interpolating  scaling  function  $  constructed  from  (f>.  We 
have 

■ah- 


Writing  9(u>)  =  Ylt ~-<x>  ^kBkuJ,  we  see  that  9(2  3u>)  is  the  Fourier  transform  of  the 
distribution 

OO 

(32)  Qj(x)  =  £  9k8(x-k 2~k). 

k—~  oo 


Thus  $  is  derived  from  (f>  via  an  iterated  convolution: 

4>M  =  ((*S.e)W)M> 

which  in  view  of  (32)  corresponds  to  an  iterated  averaging  of  pointwise  function  values. 
This  will  tend  to  have  a  smoothing  effect. 


12.  The  interpolating  multiresolutions  associated  with  the  spline  mul¬ 
tiresolutions.  In  this  section  we  will  construct  interpolating  multiresolutions  from  the 
spline  multiresolutions  of  Lemarie  and  Battle.  The  spline  multiresolution  of  order  m 
has  the  following  properties  (see  [4]): 

1.  The  scaling  function  <f>  and  the  wavelet  rjj  are  both  Cm~2  splines.  On  each 
interval  [k,k  +  1],  k  6  Z,  <f>  is  a  polynomial  of  degree  m  —  1,  while  on  each 
interval  [k/ 2,  ( k  +  l)/2],  k  €  Z,  ijj  is  a  polynomial  of  degree  m  —  1. 

2.  Both  <j)  and  are  exponentially  decaying.  There  exist  C  >  0  and  e  >  0  such 
that 


W*)l,lfWI....,|#(”-J)(*)|  <  Ce- W 
|V.WI,l'/-'(UI,....|*/’(m“2)(U|  <  Ce- H 
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Let  Hm  denote  the  scale-transition  filter  associated  with  the  Lemarie- Battle  spline 
multiresolution  of  order  m.  It  has  the  form 


Hm(u) 


^2 m(kQ 
\J  22m£2m(2cu)  ’ 


where 


.  1  d2m  2  1  _  Pm{ COS2  |) 

-  (2m-  l)\duj2m-2  4  sin2  f  ~  sin2m  f  ' 

Pm  is  a  polynomial  of  degree  m  and  Pm(c os2  ^)  >  0  for  all  u  6  R.  The  scaling  function 
(f>  and  the  wavelet  ip  associated  with  the  spline  multiresolution  of  order  m  are  given  by 

—  uj  m  (S2m(t<;))  2 

+  ir)^(^). 

For  odd  m  we  can  construct  an  interpolating  multiresolution  using  the  method 
described  in  Section  11.  Set 


a(w) 

b(u) 


Hm{u)  1 2  =  COS 


2m  W  Pm(cOS  2) 
2  Pm(cos2  lo)  ’ 


m  U)  .  m  UJ 

cosm  -  sin  - 

2 


P to  (cos2  f)PTO(sin2  f) 
P2(cos2u;) 


Then  Hi{u>)  =  a(w)  +  i6(w)  defines  an  interpolating  multiresolution  which  inherits 
the  smoothness,  decay,  and  vanishing  moment  properties  of  the  spline  multiresolution. 
Applying  Theorem  11.1,  we  obtain  the  following. 

PROPOSITION  12.1.  Suppose  Hi(u>)  is  defined  as  above.  Then  H\{ui)  defines  a 
regular  multiresolution  for  which  the  scaling  function  $  and  ivavelet  'F  have  the  following 
properties. 

1.  The  scaling  function  $  and  the  wavelet  'F  are  both  real-valued,  Cm~ 2  functions. 

2.  $  is  continuous  and  $(k)  =  6(k). 

3.  Both  $  and  *F  are  exponentially  decaying.  There  exist  C'  >  0  and  e'  >  0  such 
that 


i«wi,in»)i . 


<  C'e~e'w 

<  C'e 


4.  Both  $  and  *F  have  vanishing  moments: 


Jdxxk$(x)  —  0 
J dx  xk' F(.t)  = 


0 


for  k  =  1, ...  ,m  —  1,  and 
for  k  =  0, . . . ,  m  —  1. 
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The  scaling  function  <fm(x)  of  the  spline  multiresolution  of  order  m  approximates 
a  lowpass  filter  for  the  frequency  interval  [ — 7r,  7r],  which  is  the  Nyquist  frequency  asso¬ 
ciated  with  the  sampling  interval  Ax  =  1.  In  fact,  as  m  — >  oo,  <fm  (k'O  *  X[— 7r,7r](^)?  fhe 
characteristic  function  for  [ — 7r,  7t]  (see  [4]).  (Of  course,  the  penalty  for  such  spectral 
localization  is  that  as  m  — >  oo,  the  function  (f>(x)  is  less  localized  in  space.) 

The  interpolating  scaling  function  has  the  property  that  <f>m(o>)  =  ,  so 

it,  too,  has  most  of  its  frequency  content  concentrated  in  the  interval  [ — 7r,  7r].  This  may 
make  4>m  an  interesting  interpolant  for  sampled  functions,  where  there  is  an  underlying 
assumption  that  the  data  is  band-limited. 

Figures  1  and  2  show  plots  of  $  and  'P  for  the  case  m  =  3. 

13.  The  interpolating  multiresolutions  associated  with  the  compactly 
supported  multiresolutions.  We  can  also  build  interpolating  multiresolutions  from 
the  multiresolutions  of  Daubechies  [2],  for  which  the  associated  scaling  function  <f>  and 
wavelet  t/>  are  compactly  supported.  The  Daubechies  multiresolution  of  order  N  has 
the  following  properties: 

1.  <f>  and  are  compactly  supported. 

2.  <f>  and  0  are  smooth,  in  the  sense  that  there  exists  >  0  such  that 

(1  +  |^|  and  (1  +  M)1+ajv  |^|  6  L°°. 

The  scale-transition  filter  Hn(lj)  for  the  Daubechies  multiresolution  of  order  TV  has 
the  form 

Hn(u)=  [e'fcos  ^\N  Q(e% 
where  Q  is  a  polynomial  such  that 


.2  N~l  ( 

<2(e“)|  =  £ 

k=0  \ 


N-l+k 

k 


sin2fc  —  +  sin2  —  R  (  -  costu 


and  R  is  a  member  of  a  restricted  class  of  odd  polynomials.  For  our  construction  we 
will  choose  the  simplest  possibility,  R  =  0. 

For  odd  TV  we  start  with  to  define  an  interpolating  multiresolution.  Define 


a(u)  =  \Hn(u)\2  =  cos2A  ^  Q{e,ul)  , 


b(u)  =  cosN^s\nN  |  |Q(e*‘“)|  |g(e!(^)| . 

Now  we  can  now  appeal  to  Theorem  11.1  to  define  an  interpolating  multiresolution  from 
H\(u>)  =  a(u)  -f  ib( u>). 

PROPOSITION  13.1.  Suppose  Hi(u)  is  defined  as  above.  Then  defines  a 

regular  multiresolution  for  which  the  scaling  function  $  and  wavelet  have  the  following 
properties. 
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1.  If  a n  is  the  exponent  in  Property  2  of  the  compactly  supported  multiresolution , 
then 


(1  +  \io\)1+aN  $ 


and  (1  +  M)1+a" 


e  L°°. 


2.  $  is  continuous  and  <&(k)  =  8(k). 

3.  $  and  'll  both  have  vanishing  moments: 

j dx  xk<&(x)  =  0  for  k  =  1, . . . ,  N  —  1,  and 
Jdxxkty(x)  =  0  for  k  =  0, . . . ,  N  —  1. 

f.  If  m  is  an  integer  for  which  0  <  m  <  a pi,  then  $  and  'll  are  both  Cm.  Moreover, 
vp  and  their  derivatives  are  all  exponentially  decaying.  There  exist  C  >  0 
and  s  >  0  such  that 


|*(x)|,|*'(x)|,...,|*<“»(*)|  <  Ce-'W 
|*(x)|,|*'(x)|,...,|*,”>(x)j  <  Cf  W. 

Plots  of  $  and  for  the  case  N  —  3  are  given  in  Figures  3  and  4.  To  the  naked 
eye,  these  appear  to  be  a  little  smoother  than  the  plots  of  <f>  and  xf  for  the  compactly 
supported  multiresolution  of  order  3  given  in  [2].  This  is  what  we  might  expect  in  view 
of  the  averaging  effect  discussed  at  the  end  of  Section  1 1 . 
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Figure  1:  Interpolating  scaling  function 
built  from  spline  multiresolution,  m  =  3 


Figure  2:  Associated  wavelet  built 
from  spline  multiresolution,  m  =  3 
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Figure  3:  Interpolating  scaling  function  built 
from  Daubechies  multiresolution,  N  =  3 


Figure  4:  Associated  wavelet  built  from 
"Daubechies  multiresolution,  N  =  3 


